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ABSTRACT 

We describe an automated method for detecting clusters of galaxies in imaging 
and redshift galaxy surveys. The Adaptive Matched Filter (AMF) method utilizes 
galaxy positions, magnitudes, and — when available — photometric or spectroscopic 
redshifts to find clusters and determine their redshift and richness. The AMF can be 
applied to most types of galaxy surveys: from two-dimensional (2D) imaging surveys, 
to multi-band imaging surveys with photometric redshifts of any accuracy (2lD), to 
three-dimensional (3D) redshift surveys. The AMF can also be utilized in the selection 
of clusters in cosmological N-body simulations. The AMF identifies clusters by finding 
the peaks in a cluster likelihood map generated by convolving a galaxy survey with 
a filter based on a model of the cluster and field galaxy distributions. In tests on 
simulated 2D and 2^D data with a magnitude limit of r' ~ 23.5, clusters are detected 
with an accuracy of Az ~ 0.02 in redshift and ~10% in richness to z ^ 0.5. Detecting 
clusters at higher redshifts is possible with deeper surveys. In this paper we present 
the theory behind the AMF and describe test results on synthetic galaxy catalogs. 



1. Introduction 



Clusters of galaxies — the most massive virialized systems known — provide powerful tools 



in the study of cosmology: from tracing the large-scale structure of the universe ( Bahcall 1988 



Huchra et al 1990|, postman, Huchra & Geller 1992 


, Dalton et al 1994, 


Peacock & Dodds 1994 


and references therein) to determining the amount of dark matter on Mpc scales (Zwicky 1957 


Tyson, Valdes & Wenk 199(1 


Kaiser & Squires 1993, 


Peebles 1993, 


Bahcall, Lubin &: Dorman 



1995, Carlberg et al 1996| ) to studying the evolution of cluster abundance and its cosmological 



implications (Evrard "19891 , Peebles, Daly fc Juszkiewicz 1989| , Penry et alT99^ , [Eke, Cole &] 
Frenk 19961 , pahcall, Fan k Cen 199~7| , |Carlberg et al 199% |Oukbir fc Blanchard 19971 ). The above 
studies place some of the strongest constraints yet on cosmological parameters, including the 
mass-density parameter of the universe, the amplitude of mass fluctuations at a scale of 8 h _1 Mpc 
and the baryon fraction. 
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The availability of complete and accurate cluster catalogs needed for such studies is limited. 
One of the most used catalogs, the Abell catalog of rich clusters ( [Abell (1958) , and its southern 



counterpart Abell, Corwin & Olowin 1989), has been extremely useful over the past four decades. 
This catalog, which contains ~4000 rich clusters to z ^ 0.2 over the entire high latitude sky, 
with estimated redshifts and richnesses for all clusters, was constructed by visual selection from 
the Palomar Sky Survey plates, using well-defined selection criteria. The Zwicky cluster catalog 
( Zwicky et al 1961-68| ) was similarly constructed by visual inspection. 



The need for new, objective, and accurate large-area cluster catalogs to various depths is 
growing, following the important use of clusters in cosmology. Large area sky surveys using CCD 
imaging in one or several colors, as well as redshift surveys, are currently planned or underway, 
including, among others, the Sloan Digital Sky Survey (SDSS). Such surveys will provide the 
data needed for constructing accurate cluster catalogs that will be selected in an objective 
and automated manner. In order to identify clusters in the new galaxy surveys, a robust and 
automated cluster selection algorithm is needed. We propose such a method here. 

Cluster identification algorithms have typically been targeted at specific surveys, and new 



algorithms have been created as each survey is completed. |Abell (1958] was the first to develop a 
well-defined method for cluster selection, even though the identification was carried out by visual 
inspection (see, e.g., McGill fc Couchman 199~C| for a analysis of this method). Other algorithms 



have been created for the APM survey ( palton et al 1994 |Dalton, Maddox fe Sutherland 1997] 



see |Schuecker fc Bohringer 1998 for a variant of this method), the Edinburgh-Durham survey 



(ED; [Lumsden et al 1992j ), and the Palomar Distant Cluster Survey (PDCS; postman et al 1996| ; 
see also Kawasaki et al 1997 for a variant of this method; and |Kleyna et al 199£ for an application 



this method to finding dwarf spheroidals). All the above methods were designed for and applied 
to two-dimensional imaging surveys. 

In this paper we present a well defined, quantitative method, based on a matched filter 
technique that expands on some of the previous methods and provides a general algorithm that can 
be used to identify clusters in any type of survey. It can be applied to 2D imaging surveys, 2^D 
surveys (multi-band imaging with photometric redshift estimates of any accuracy), 3D redshift 
surveys, as well as combinations of the above (i.e. some galaxies with photometric redshifts and 
some with spectral redshifts). In addition, this Adaptive Matched Filter (AMF) method can be 
applied to identify clusters in cosmological simulations. 

The AMF identifies clusters by finding the peaks in a cluster likelihood map generated by 
convolving a galaxy survey (2D, 2^D or 3D) with a filter which models the cluster and field galaxy 
distribution. The peaks in the likelihood map correspond to locations where the match between 
the survey and the filter is maximized. In addition, the location and value of each peak also gives 
the best fit redshift and richness for each cluster. The filter is composed of several sub-filters that 
select different components of the survey: a surface density profile acting on the position data, a 
luminosity profile acting on the apparent magnitudes, and, in the 2^D and 3D redshift 
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cut acting on the estimated redshifts. 

The AMF is adaptive in three ways. First, the AMF adapts to the errors in the observed 
redshifts (from no redshift information (2D), to approximate (2^D) or measured redshifts (3D)). 
Second, the AMF uses the location of the galaxies as a "naturally" adaptive grid to ensure 
sufficient spatial resolution at even the highest redshifts. Third, the AMF uses a two step approach 
that first applies a coarse filter to find the clusters and then a fine filter to provide more precise 
estimates of the redshift and richness of each cluster. 

We describe the theory of the AMF in §2 and its implementation in §3. We generate a 
synthetic galaxy catalog to test the AMF in §4 and present the results in §5. We summarize our 
conclusions in §6. 

2. Derivation of the Adaptive Matched Filter 

The idea behind the AMF is the matching of the data with a filter based on a model of 
the distribution of galaxies. The model describes the distribution in surface density, apparent 
magnitude, and redshift of cluster and field galaxies. Convolving the data with the filter produces 
a cluster probability map whose peaks correspond to the location of the clusters. Here we describe 
the theory behind the AMF. We present the underlying model in §2.1, the concept of the cluster 
overdensity in §2.2, two likelihood functions derived under different assumptions about the galaxy 
distribution in §2.3, and the extension of the likelihood functions to include estimated redshifts in 
§2.4. 

2.1. Cluster and Field Model 

The foundation of the AMF is the model of the total number density of galaxies per solid 
angle (dQ = 2ir9d9) per apparent luminosity (dl) around a cluster at redshift z c with richness A 

n mo dei(#, l\ z c )d£ldl = [rif(l) + An c (9, 1; z c )]dQdl , (1) 

where 9 is the angular distance from the center of the cluster, and n/ and An c are the number 
densities due to the field and the cluster. The field number density is taken directly from the global 
"number-magnitude" per square degree relation. The cluster number density is the product of a 
projected density profile (see Appendix A) and a luminosity profile, both shifted to the redshift of 
the cluster and transformed from physical radius and absolute luminosity to angular radius and 
apparent luminosity (i.e. flux) 

An c (9,l;z c ) = AS c (r) (^) (§) (2) 
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where the projected radius r and absolute luminosity L at the cluster redshift are given by 

0d c 



r(6,z c 



l + z, 
i2 



L(l,z c ,n) = AirD\z c ,n)l . (3) 

Here d c is the angular size distance corresponding to z c , and D is the luminosity distance (with 
K correction) of a galaxy of spectral type \x (e.g. elliptical, spiral or irregular) (see Appendix B). 
The factor of 1 + z c in the radius relation converts from comoving to physical units (the cluster 
profile can be defined with either comoving or physical units). For the cluster luminosity profile a 
Schechter function is adopted 

</> c (L)dL oc (L/L*)- a e~ L/L *d(L/L*), (4) 
where a ~ 1.1 ( [Schechter 1976| , Binggeli, Sandage fc Tammann 1988| , Loveday et al 199^ ). 



The model is completed by choosing a normalization for the radial profile S c and the 
luminosity profile cp c . The choice of normalization is arbitrary, but has the effect of determining 
the units of the richness A. We choose to normalize so that the total luminosity of the cluster is 
equal to 

L c (< r max ) = AL* 

= ^ J L(l,z c )n c (6,l;z c )dftdl 

roo 



A / S c (r)2vrrdr / (/) c (L)LdL . (5) 
Jo Jo 



max- 



The richness parameter A thus describes the total cluster luminosity (in units of L*) within r n 
For r max — 1 h ^Mpc, the richness A(< 1 h ^Mpc) relates to Abell's richness A^4 (within 1.5 
h -1 Mpc) as Na ~ §A. For example, Abell richness class > 1 clusters (Na > 50) correspond to 
A > 75 ( Bahcall &: Cen 1993| ) . Multiplying by L in equation (5) allows (f> c to be integrated down 



to zero luminosity, thus insuring that the total luminosity is equal to AL* regardless of the the 
apparent luminosity limit of the survey Z mm . The above constraint can be implemented by simply 
multiplying S c and (f) c by any appropriate constants (e.g., normalizing the first integral to one and 
setting the second integral equal to L*). 



2.2. Cluster Overdensity 

Clusters, by definition, are density enhancements above the field. To quantify this we 
introduce the cluster overdensity AA, which is the sum of the individual overdensities of the 
galaxies A<5j. A and 5i are defined subsequently. 

For a given cluster location on the sky let 0^ and Z, be the angular separation from the cluster 
center and the apparent luminosity of the ith galaxy, respectively. For a specific cluster redshift z c 
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we only need consider galaxies inside the maximum selected cluster radius (Appendix A) 

n < y \_ (1 + z c )r m3jX 

The apparent local overdensity at the location of the ith. galaxy is given by A<5j, where 



(6) 



< = n c (6i,li;z c )dndl 
1 n) n f (k)dndl ' { ' 

The apparent overdensity of the cluster as measured from the data is 

AA data = A]>>, (8) 

i 

where the sum is carried out over all galaxies within mSLX . The cluster overdensity can also be 
calculated from the model 

POO ^max 

AA mode i = A/ / 6(6,l;z c )An c (6,l;z c )dndl . (9) 
Equating Adata with A mo dd it is possible to solve for A 

^ = fr L - (io) 

J n c 

The term in the denominator is simply the model cluster overdensity when A = 1. As we shall see 
in the subsequent sections, the positions of clusters, both in angle and in redshift, correspond to 
the locations of maxima in the cluster overdensity. 



2.3. Likelihood Functions 

Having defined the model we now discuss how to find clusters in a galaxy catalog and 
determine the best fit redshift and richness z c and A for each cluster. The basic scheme is to 
define a function C, which is the likelihood that a cluster exists, and to find the parameters that 
maximize £ at a given position over the range of possible values of z c and A. This procedure is 
carried out over the entire survey area and generates a likelihood map. The clusters are found by 
locating the peaks in the likelihood map. The map grid can be chosen by various means, such as 
a uniform grid or by using the galaxy positions themselves. 

A variety of likelihood functions can be derived, depending on the assumptions that are made 
about the distribution of the galaxies. The AMF uses two likelihood functions whose derivations 
are given in Appendix C. A cluster is identified, and its redshift and richness, z c and A, are found 
by maximizing C. Typically this is accomplished by first taking the derivative of C with respect 
to A and setting the result equal to zero 

*£-<>■ <"> 
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One can compute A from the above equation either directly or numerically. The resulting value of 
A can be inserted back into the expression for C to obtain a value for C at the specified redshift. 
Repeating this procedure for different values of z c it is possible to find the maximum likelihood 
and the associated best value of the cluster redshift, as well as the best cluster richness at this 
redshift. 

The first likelihood, which we call the coarse likelihood, assumes that if we bin the galaxy 
catalog, then there will be enough galaxies in each bin that the distribution can be approximated 
as a Gaussian. This assumption, while not accurate, provides a coarse likelihood function £ C oarse 
which is linear, easy to compute and corresponds to the apparent overdensity 

■^coarse — -^-coarse ^ ' $i i (1^) 

i 

where 

Acoarsc = f5(e,l)n c (e,l)cMdl (13) 

(see Appendix C). 

The second likelihood, which we call the fine likelihood, assumes that the galaxy count in a 
bin can be modeled with a Poisson distribution — an assumption which is nearly always correct. 
The resulting likelihood function £fi ne is non-linear and requires more computations to evaluate, 
but should provide a more accurate estimate of the cluster redshift and richness. The formula for 

•£finc IS 

£ fine = -A fine iV c + M 1 + Afine<y , (14) 

i 

where Ag ne is computed by solving 



N, 



c 



+ (15) 



and N c is the number of galaxies one would expect to see in a cluster with A = 1 (see Appendix 
C). 



2.4. Including Estimated Redshifts 

So far we have considered cluster selection in purely 2D imaging surveys with no estimated 
redshift information. The inclusion of estimated galaxy redshifts, as in 2^D or 3D surveys, 
should improve the accuracy of the resulting cluster redshifts. In addition, estimated redshifts 
can separate background galaxies (noise) from cluster galaxies (signal) more easily, allowing the 
detection of considerably poorer clusters than if estimated redshifts were not used. 

Galaxy redshift information can be obtained, for example, from multi-color photometric 
surveys (via photometric redshifts estimates) or from direct spectroscopic determination of galaxy 
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redshifts. We now discuss how to extend the AMF described above to include such redshifts. The 
galaxy redshifts that we use within a given survey can range from very precise measured redshifts, 
to only approximate photometric redshifts (with varied accuracy), to no redshift information — all 
in the same analysis (i.e. adapting to different redshift information for each galaxy in the survey). 
The available redshifts provide essentially a third filter, in addition to the spatial and luminosity 
filters used in the 2D case. In practice, we use the estimated redshift information of each galaxy 
to narrow the window of the AMF search following the same basic method described above. 

One of the benefits of laying down the theoretical framework for the AMF is the easy means 
by which estimated redshifts can now be included. Let Zj and a\ be the estimated redshift and 
redshift error of the ith galaxy, and let the factor w determine how wide a region we select around 
each z c for the likelihood analysis. Inclusion of the additional redshift information is accomplished 
by simply limiting the sums in £ coa rse and £fi ne to those galaxies that also satisfy \zi — z c \ < wa l z . 
Note that this procedure allows the usage of redshift information with variable accuracy in the 
same survey — i.e., some galaxies with measured redshifts, some with estimated redshifts and some 
with no redshifts at all. 

The richness A coarse is dominated by the cluster galaxies; as long as w is sufficiently large (e.g. 
w ~ 3) then the resulting richnesses will be unbiased. However, if the redshift errors are large 
(i.e., a large fraction of the depth of the galaxy catalog), it may be desirable to use a smaller value 
of w « 1. In this case a small correction needs to be applied to the richness to account for the 
small fraction of cluster galaxies that are eliminated by the redshift cut. If the redshift errors are 
Gaussian, the desired correction can be obtained by multiplying the predicted cluster overdensity 
by the standard error function 

/2 r i,2 , 

\ - / e"5* dt . (16) 

V 7T Jo 

A similar issue arises with Afi ne when the field galaxies are eliminated by using estimated 
redshifts. To obtain the correct value of Afi ne requires modeling the contribution of the eliminated 
field galaxies to equation (15), which can be done either analytically or numerically via Monte 
Carlo methods. 



3. Implementation 

The likelihood functions derived in the previous section represent the core of the AMF cluster 
detection scheme. Both likelihood functions begin with picking a grid of locations over the survey 
area for which C and A are computed over a range of redshifts. Figures [I] and || show the functions 
C(z) and A(z) for a position on the sky located at the center of a A = 100 cluster (~ | the richness 
or luminosity of Coma) at z = 0.35 (the details of the test data are given in the next section). 
Figure [l] shows the results when no redshift estimate exists (i.e., a z — ► oo), and figure § shows 
what happens when using estimated redshifts with 0.03 < o z < 0.06. 
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Figures || and |2| illustrate three of the basic properties of these likelihood functions, which 
were discussed in the previous section: 1. the dramatic effect of including estimated redshifts on 
sharpening the peak in C, which makes finding clusters much easier; 2. the qualitative difference 
in the form of Afi ne (z) and A coarse (z) in figure [2] that arises from the fact that A coarse is a simple 
function of £ C oarse while Afi ne is a complicated function of £fi ne ; 3. the shortward bias in the 
cluster redshift as computed from £ C oarse in 2D (figure ||) , which is a general bias in the 2D coarse 



likelihood function (for a discussion of how to correct for this bias see [Postman et al 1996|) . In 
the 3D case, the coarse likelihood of rich clusters has comparable accuracy to the fine likelihood. 
For poorer clusters (A < 50), the coarse likelihood yields higher richness estimates than the true 
values; this is a result of the Gaussian assumption used for deriving the coarse likelihood which 
worsens for poorer clusters. 

Implementation of the AMF cluster selection consists of five steps: (1) reading the galaxy 
catalog, (2) defining and verifying the model, (3) computing £ coa rse over the entire galaxy survey 
over a range of redshifts, (4) finding clusters by identifying peaks in the £ C oarse map, and (5) 
evaluating £fi ne and obtaining a more precise determination of each cluster's redshift and richness. 



3.1. Reading the Galaxy Catalog 

The first step of the AMF is reading the galaxy catalog over a specified region of the sky. 
The galaxy catalog consists of four to six quantities for each galaxy, i: the position on the sky 
RAj and DECj, the apparent luminosity (i.e. flux) in specific band li, the Hubble type used for 
determining the K correction (e.g., E, Sa or Sc), the estimated redshift Z{ and estimated error 
a\. In the case of a single band survey, where no photometric redshifts are possible, Z{ and a\ 
will not exist. In addition to these local quantities, the following global quantities are computed 
from the catalog: the apparent luminosity limit ^ m i n , the mean estimated error a z . Finally, it is 
necessary to set the minimum and maximum cluster redshifts for the cluster search z m { n and z mSjX 



3.2. Model Definition and Verification 

The primary model components which need to be specified are: the cluster radial profile 
£ c (r), the cluster luminosity profile (f) c (L), a normalization convention for the cluster that sets 
the units of richness A, and the field number density versus apparent luminosity distribution 
nf(l). In addition, a specific cosmological model needs to be chosen along with K corrections for 
each Hubble type, from which the angular distance d(z) and luminosity distance D(z, fi) can be 
computed. The observed galaxy catalog needs to be tested for its consistency with the model 
parameters (such as nf(l)). While this is unnecessary for the simulated galaxy catalog used below 
(where the model used to find clusters is the same as the one used to generate the catalog), we 
show some of the consistency checks for illustration purposes (see Appendix D). 
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3.3. £ C oarse Mapping 

Mapping out the coarse likelihood function begins with picking a grid covering the survey area. 
The most straightforward choice is a uniform grid over the area covered, which is conceptually 
simple and makes finding peaks in the likelihood easy. However, a uniform grid needs to be 
exceedingly fine in order to ensure adequate resolution at the highest redshifts, and leads to 
unacceptable computer memory requirements. Another choice for the grid locations is to use the 
positions of the galaxies themselves. The galaxy positions are a "naturally" adaptive grid that 
guarantees sufficient resolution at any redshift while also eliminating unnecessary points in sparse 
regions. For this reason we choose to use the galaxy positions as the grid locations. 

At each grid location (i.e. each galaxy position) £ coa rsc is evaluated over a range of redshifts 
from z m \ a to z max . The number of redshift points is set so there is adequate coverage for the given 
value of a z . Finding the maximum of the likelihood sets the values of the likelihood, redshift and 
richness at this point: -Ccoarse' z coarse> an d -^-coarse- When this process has been completed for all 
the galaxies, the result is an irregularly gridded map £* oarse (RAj, DECj, z* oarse ). The peaks in this 
map correspond to the locations of the clusters. 



3.4. £ C oarsc Cluster Selection 

Finding peaks in 3D regularly gridded data is straightforward. Finding the peaks in the 
irregularly gridded function £* oarse (RAj, DECj, z* oarse ) is more difficult. There are several possible 
approaches. We present a simple method which is sufficient for selecting rich clusters. More 
sophisticated methods will be necessary in order to find poor clusters that are close to rich clusters. 

As a first step we eliminate all low likelihood points £ coarse < £ C ut> where £ cut is the nominal 
detection limit, which is independent of richness or redshift. £ cut can be estimated from the 
distribution of the hoarse values (figure ||). The peak in the distribution corresponds to the 
background while the long tail corresponds to the clusters. Assuming the background is a Gaussian 
whose mean can be estimated from the location of the peak and whose standard deviation is 0.43 
times the Full- Width-Half-Maximum (FWHM), then a given value of £ cu t will lie N a standard 
deviations from the peak 

a ~ 0.43 FWHM ' 1 ' 

The values of C cu t shown in figure |3] were chosen so that N a ~ 5, which was sufficiently high that 
no false detections occurred. To first order, £ cut oc a z /z, where z is the average depth of the 
survey and is a function of l m i n - 

Step two consists of finding the largest value of £* oarse , which is by definition the first and 
most overdense cluster. The third step is to eliminate all points that are within a certain radius 
and redshift of the cluster. Repeating steps two and three until there are no points left results 
in a complete cluster identification (above £ C ut); with a position, redshift and richness (oc total 
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luminosity) estimate for each cluster. 



3.5. £fi nc Evaluation 

The angular position, redshift and richness of the clusters determined by the £ C oarse selection 
are adequate, but can be complemented by determining the redshift and richness from £fi nc . 
Recall that £fi nc requires 10 to 100 times as many operations to evaluate as compared with £ C oarse> 
and applying it to every single galaxy position can be prohibitive. Evaluating £fi nc on just the 
clusters found with £ C oarse is trivial and worth doing as it should provide more accurate estimates 
of the cluster redshift and richness because of the better underlying assumptions that went into 
its derivation. Thus, the final step in our AMF implementation is to compute Zg nc and Ag ne using 
ne on each of the cluster positions obtained with £ coarse . 



3.6. Implementation Summary 

We summarize the above implementation in the following step-by-step list. 

1. Read galaxy catalog 

• Read in RAj, DECj, k, Zi, a\ and fii for each galaxy. 

• Pick Z m i n (survey limit), ^mm 

and z max (cluster search limits), and compute average 

distance error a z . 

2. Model definition 

• Choose S c , 4> c , rif, d(z c ) and D(z,/j,). 

• Normalize S c and <p c . 

• Verify galaxy distributions in I and z with those predicted by model. 

3. ^coarse mapping 

• For each galaxy location (RAj, DEQ) choose a range of redshifts within z m i n and z max 
with a step size no larger that one half u z . 

• Compute £ coa rsc (eq. 12) and A coarse (eq. 13) for each redshift. Set £* oarse , A* oarsc and 
hoarse equal to the values at the maximum of £ C oarse • 

4. Find peaks in the £ C oarsc map 

• Compute C cu t from the C l coaise distribution (e.g. 5-a cut). 

• Find all local maxima in £* oarse (RAj, DECj, z* oarse ) where £ l coarse > £ cut ; these are the 
clusters. 
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5. Refine cluster redshift and richness estimates with £fi ne 

• At the RA and DEC of each cluster found with £ coarse evaluate £fi ne over the same 
range of redshifts within z m \ n and z max . 

• Compute £fi ne (eq. 14) and Afi ne (eq. 15) for each redshift. Set £g ne , Ag ne and z\ ne 
equal to the values at the maximum of £fi ne . These provide the best estimates for the 
cluster richness and redshift. 

In the next section we describe a simulated galaxy catalog used to test the above AMF 
implement ation . 



4. Simulated Test Data 

Our test data consists of 72 simulated clusters with different richnesses and redshifts placed in 
a simulated field of randomly distributed galaxies (for a non-random distribution of field galaxies 
see section 5). The clusters range from groups to rich clusters with total luminosities from 10 
L* to 300 L* (corresponding to Abell richness counts of approximately 7 to 200, or richness 
classes <C0 to ~4) and are distributed in redshift from 0.1 to 0.5. The luminosity profile is a 
Schechter function with a = 1.1. The radial profile is a Plummer law given in Appendix A with 
''max = 1 h ^Mpc and r core = 0.1 r max . The number and absolute luminosity of the field galaxies 
were generated from a Schechter function, using the field normalization (j)* = 1.08 x 10 -2 h 3 Mpc~ 3 



flLoveday et al 1992|) 



Three different Hubble types were used, E, Sa and Sc, with K corrections taken from Poggianti 



1997 , Each galaxy in a cluster was randomly assigned a Hubble type so that 60% were E, 30% 
were Sa, and 10% were Sc; each galaxy in the field is randomly assigned a Hubble type so that 
40% were E, 30% were Sa, and 30% were Sc. Knowing the redshift of each galaxy, its absolute 
luminosity and Hubble type, we can compute the apparent luminosity. All galaxies with apparent 
luminosity below r' ~ 23.5 (the anticipated limit of the SDSS) were eliminated. This limit resulted 
in a field number density of ~5000 galaxies per square degree. 

To facilitate the subsequent analysis and interpretation of the results, the clusters were placed 
on an 8 by 9 grid. The cluster centers were separated by 0.4 degrees, which results in the test data 
having dimensions of 3.2 by 3.6 degrees. The distribution of the cluster galaxies in RA and DEC 
is shown in figure where each column has the same richness while each row of clusters is at the 
same redshift (see figure ||). From left to right the richnesses are A = 10, 20, 30, 40, 50, 100, 200, 
and 300. From bottom to top the redshifts are z c = 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, and 
0.5. In all, the clusters contained some 30,000 galaxies, over half of which lie in the richness 200 
and 300 clusters. The randomly generated field contained approximately 50,000 galaxies. The 
distribution of all the galaxies in RA and DEC is shown in figure ||. 

So far, the redshifts for the galaxies are exact. If we assume that redshift errors are Gaussian, 
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then we can easily generate offsets to the true positions if we know the standard deviation of 
the distribution o z . Not all the galaxies will have the same estimated redshift error a l z , and, in 
the case of photometric redshifts, these values can be expected to vary by about a factor of two 
( Connolly et al 1995| ; Yee 1998). We model this effect by first randomly generating the estimated 
redshift errors from a uniform distribution over a specified range (e.g., 0.03 < a\ < 0.06). The 
offsets from the true redshifts are then randomly generated from Gaussian distributions with 
standard deviations corresponding to o\ values. The estimated redshifts are computed by adding 
the offsets to the true redshift. Thus, a dataset with o z ~ 0.05 refers to 0.03 < a\ < 0.06. 



5. Results and Discussion 

To test our AMF implementation we ran it on the above simulated galaxy catalog for three 
different error regimes: a z ~ 0.05, a z ~ 0.15 and a z — > oo. A contour plot computed from the 
maximum likelihood points -Qoarse is shown in figure [7j which indicates that the coarse filter does 
a good job of finding the angular positions of the clusters with no false detections. The resulting 
cluster centers have an accuracy that is within one core radius of the true cluster center. The 
redshifts and richnesses obtained by applying £g ne to the clusters found with £ coa rse are shown in 
figures ||, ||, and [U]. The boxes denote the true redshift and richness values of the input clusters. 
The short lines connect the inputs with the outputs (i.e. the values obtained with £fi no ). Smaller 
lines indicate more accurate redshift and richness determinations. The long curve indicates the 
expected detection limit for the value of £ cut used. As expected, the number of clusters detected 
and their accuracy decrease as a z increases. However, nearly all the clusters with A ^ 100 
(corresponding roughly to Abell richness class ^ 1) are detected out to redshifts of 0.5 which is 
the effective limit of the survey. In a deeper survey it will be possible to detect clusters at higher 
redshifts. 

The errors in the redshift and richness estimates of all detected clusters is presented in 



figures 11 and 12. A summary of these results is shown in table 1. In general, the AMF finds 
clusters with an accuracy of Az ~ 0.02 in redshift and ~10% in richness. Including distance 
information lowers the background and results in a substantial improvement in the detection of 
poorer clusters. Thus, many more clusters are detected when a z ~ 0.05 as compared to a z —* oo. 
These additional clusters are all poorer and thus have higher errors, which explains why the 
average errors do not significantly change with a z . 

Six additional tests were conducted on the a z ~ 0.05 case to check the robustness of the 
results. The first test explored the effect of spatial structure in the background distribution 
of galaxies by using positions taken from an N-body simulation (instead of using randomly 
distributed field galaxies). The second test investigated the effect of using different K-corrections. 
The next four tests explored the effect of changing different parameters of the cluster model: a in 
the Schechter luminosity function, n in the Plummer cluster density profile, the core radius (r core ) 
and the maximum radius of the cluster profile (r max ). The results of all these tests are summarized 
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Table 1: Summary of AMF tests on simulated data 



Input 


Likelihood 


Output 


Output 


error 


cutoff 


z error 


A error 




■^cut 


a(Az) 


er(AA/A) 


0.05 


40 


0.014 


0.13 


0.15 


100 


0.025 


0.12 


oo 


300 


0.023 


0.14 



in Table 2. 



The N-body simulation contained 128 3 dark matter particles in a 200 h _1 Mpc box (|Xu 1995|) 



with sufficient spatial resolution to resolve cluster cores. The final z = ouput of the simulation 
was "stacked" in a non-repeating fashion ( pott 1997 ) to create a simulated field out to z = 0.6. 



Each particle was then assigned a luminosity in the same manner as described earlier for the 
uniform background. The 72 test clusters were placed in the N-body background. The coarse 
filter was run on these data using the same parameters as in the uniform case. All the clusters 
detected with the uniform background were also detected with the N-body background. Next, 
the fine filter was run on these clusters. There was no change in the estimated redshifts of the 
clusters. The richness estimates showed a slight increase (~ 20%) in their variance, which is due 
to the fluctuations in the background density. 

In a real survey, it is unlikely that all the galaxies will have correctly assigned Hubble types. 
An error in the Hubble type primarily affects the K-correction. To test the effects of incorrect 
Hubble types we ran the fine filter with an assumption that all the galaxies were Ellipticals (E) 
and then again assuming that all the galaxies were Spirals (Sc). These changes produced no 
significant change in the estimated redshift or richness of the clusters out to z ~ 0.5. 

In the real Universe, clusters can not be described by a single set of model parameters. Two 
tests of each of four model parameters were carried out to look at the errors produced by using 
an AMF cluster model with parameters different from the clusters in the data. In each case 
the fine filter was run on data with a z ~ 0.05 and the differences in the estimated redshift and 
richness were examined. The results are summarized in Table 2. None of the changes in any of the 
parameters significantly affected the estimated redshifts of the clusters. The largest effect on the 
richness A occurred with changing the parameter a in the Schechter luminosity function of the 
cluster. This induced a bias in A at small redshift; the bias decreases with redshift because the 
galaxies at high redshift are not part of the faint-end Schechter luminosity function. Changing the 
slope of the cluster density profile had only a small effect on A. Changing the core radius had no 
significant effect on A. As expected, A increases with increasing r max . 

The CPU and memory requirements of the AMF are dominated by the £ C oarse evaluation. 
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Table 2: Robustness tests on simulated data with a z ~ 0.05. A "no change" entry 
means that any difference in the estimated redshift or richness was within the nominal 
errors quoted in Table 1 (i.e., Az ~ 0.02 and AA/A ~ 0.1). Biases are given relative 
to the nominal value (e.g., a bias of 1.1 A implies that the new estimated richnesss is 
1.1 times the nominal estimated richness). All biases are independent of redshift and 
richness unless stated otherwise. 



Model Parameter 


New 




Effect on 


Effect on 


(nominal value) 


Value 




Estimated z 


Estimated A 


Background distribution 


N-body 


no change 


slight (~ 20%) increase in variance 


(uniform) 










K correction 


all E 




no change 


no change 


(E, Sa and Sc) 


all Sc 




no change 


no change 


Cluster lum func slope 


-0.8 




no change 


1.5A bias at z ~ 0.1, 


(a = -1.1) 








decreasing to no bias at z ~ 0.35 




-1.3 




no change 


0.5A bias at z ~ 0.1, 
increasing to 1.2 A at z ~ 0.5 


Cluster profile slope 


1.5 




no change 


1.1 A bias independent of z 


(n = 2.0) 


2.5 




no change 


0.95A bias independent of z 


Cluster core radius 


0.05 h- 


x Mpc 


no change 


no change 


(r corc = 0.1 h^Mpc) 


0.20 h" 


x Mpc 


no change 


no change 


Cluster max radius 


0.75 h~ 


x Mpc 


no change 


0.75 A bias independent of z 


(r max = 1.0 h^Mpc) 


1.25 h" 


x Mpc 


no change 


1.25 A bias independent of z 



The AMF required around 100 MBytes of memory and took from a few minutes (a z ~ 0.05) to a 
little under two hours (a z -► oo) using one CPU (198Mhz MIPS R10000) of an SGI Origin200. 
For example, the SDSS will be composed of approximately 1000 fields similar in size to our test 
catalog. Since finding clusters in one field is independent of all the others, it is simple to run the 
AMF on a massively parallel computer; it will be possible to run the AMF on the entire SDSS 
catalog in one day. 

6. Summary and Conclusions 

We have presented the Adaptive Matched Filter method for the automatic selection of clusters 
of galaxies in a wide variety of galaxy catalogs. The AMF can find clusters in most types of 
galaxy surveys: from two-dimensional (2D) imaging surveys, to multi-band imaging surveys with 
photometric redshifts of any accuracy (2^D), to three-dimensional (3D) redshift surveys. The 
method can also be utilized in the selection of clusters in cosmological N-body simulations. The 
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AMF is based on matching the galaxy catalog with a cluster filter that models the overall galaxy 
distribution. The model describes the surface density, apparent magnitude, and redshift of cluster 
and field galaxies. Convolving the data with the filter produces a cluster probability map whose 
peaks correspond to the location of the clusters. The probability peaks also yield the best fit 
redshift and richness of each cluster. 

The heart of the AMF is the apparent overdensity <5j which is evaluated at each galaxy position 
and has a higher value for galaxies in clusters than galaxies in the field. The apparent overdensity 
distills the entire description of the galaxy catalog into a single function. Two likelihood functions 
are derived, >C coarse and /2fi nc , using different underlying model assumptions. The theoretical 
framework of the AMF allows estimated redshifts to be included via a simple redshift filter, which 
effectively limits the sums in C coa _ rse and £fi ne to those galaxies within a window around z c . The 
maxima in the likelihood functions are used to identify cluster positions as well as their redshifts 
and richnesses. 

The AMF is adaptive in three ways. First, it adapts to the errors in the estimated redshifts. 
Second, it uses the locations of the galaxies as "naturally" adaptive grid to ensure sufficient 
resolution at even the highest redshifts. Third, it uses a two step approach that applies a coarse 
filter to initially find the clusters and a fine filter to more precisely estimate the redshift and 
richness of each cluster. 

We tested the AMF on a set of simulated clusters with different richnesses and redshifts — 
ranging from groups to rich clusters at redshifts 0.1 to 0.5; the clusters were placed in a simulated 
field of randomly distributed galaxies as well as in a non-random distribution produced by N-body 
cosmological simulations. We find that the AMF detects clusters with an accuracy of Az ~0.02 in 
redshift and ~10% in richness to z ^ 0.5 (for a simulated galaxy survey to r 1 « 23.5). In addition, 
robustness tests provide a strong indication that the AMF will perform well on observational data 
sets. Detecting clusters at even higher redshifts will be possible in deeper surveys. 

We wish to thank Guohong Xu for providing the N-body simulations used in this paper. In 
addition David Weinberg, Andy Connolly, Marc Postman, Lori Lubin, and Chris Finger provided 
helpful suggestions. We would also like to thank Dr. Micheal Kurtz for his helpful comments. 
This work was supported in part by the DoE Computational Science Fellowship Program, the 
Princeton Observatory Advisory Council, and NSF grants AST 93-15368 and GER 93-54937. 



A. 



Cluster Profile 



The projected cluster profile is derived from spherical modified Plummer law profile: 




r < r, 



r > r, 



max 



max 



(Al) 



where n«2 and the constant is used to normalize the profile. The projected profile is computed 
by integrating along the line of sight: 



which for the modified plummer profile gives 



r <r r 



E c (r) = { (i+r 2 A' c 2 orc )^ (i+rg mx Ac 2 o ro )^ (A3) 

, v > r max 



where £° is a constants set by the normalization. 



B. Luminosity Distance 

The transformation between absolute and apparent luminosities in band Ao 

L Xo =4irD\z;fi)l Xo , (Bl) 

where \x is the galaxy type (e.g. elliptical, spiral, irregular, ...) and D is the luminosity distance 
at a redshift z and is related to the angular size distance d by 

where K Xo and E Xo are the"K" and evolutionary corrections. For no evolution models E Xo = 1. 
The angular-size distance for Q + S7a = 1 cosmologies is given by 

c f 2 dz' 

d{z) = jrJ (i- fi+ fi(i + ^)3)i/2 ( B3 ) 



C. Derivation of the Likelihood Functions 

Various likelihood functions can be derived. The differences are due to the additional 
assumptions that are made about the distribution of the observations. This appendix gives the 
mathematical derivation of the two likelihood functions used in the AMF: £ coa rsc and C^ ne . Both 
derivations are conceptually based on virtually binning the data, but make different assumptions 
about the distribution of galaxies in the virtual bins. 

Imagine dividing up the angular and apparent luminosity domain around a cluster into bins. 
We assign to each bin a unique index j. The expected number of galaxies in bin j is denoted 
n moder ^ ne num ber of galaxies actually found in bin j is n^ ata . In general, the probability of 
finding n^ ata galaxies in cell j is given by a Poisson distribution 
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The likelihood of the data given the model is computed from the sum of the logs of the individual 
probabilities 

£ = ]TmP r (C2) 



C.l. Coarse Grained C 

If the virtual bins are made big enough that there are many galaxies in each bin, then the 
probability distribution can be approximated by a Gaussian 



Pj= . = =exp<^- V data ."models _ (C3) 



("data ~ n l 

27rn modcl ^ 2n Ldel 

Furthermore, if the field contribution is approximately uniform and large enough to dominate the 
noise then 



p = J _l^data modcV I ^ 

where n J model = n J j + An J c . Summing the logs of these probabilities results in the following 
expression for the coarse likelihood 

-^coarse = ^ ] hi Pj 
3 

= -1 y l n27r 4 - 1 y (n ' data ~ nJmodel)2 (C5) 

The first term does not affect the value of A which minimizes i3 CO arsc and can be dropped. In 
addition, if the bins can also be made sufficiently small, then the sum over all the bins can be 
replaced by an integral 

r _ f («data(6 1 ,/) - n mo dcl(0,O) 



--coarse 



= - / ^a(^/)~ n dcl (M)) dQdl (C6) 
J n f (0,l) 



where n<ia.ta.(9 , 1) is a sum of Dirac delta functions corresponding to the locations of the galaxies. 
Expanding the squared term and replacing ro mo dei with nj + An c yields 

r [ "data - 2 "data"/ - 2ra data An c + n) + 2n f An c + A 2 ™ 2 

^coarse = ~ / " dUdl . (C7) 

J n f 

The above expression can be simplified by setting 5 = n c /rif, dropping all expressions that are 
independent of A, and noting that A / n c is small compared to the other terms, which leaves 

£ coarsc = 2A y Si ~ A 2 J 5(0, l)n c (9, l)dKldl (C8) 
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Taking the derivative of £ C oarse with respect to A and setting the result equal to zero, we can solve 
for A coarse directly 



Acoarse J S(9 , l)n c (d , l)dQdl' 

where the denominator terms of (fr c , E c and nj is 

,2 



(C9) 



8(6,l)n c (e,l)dndl 



n: 



dVLdl 



T, 2 c (r')2nr'dr' 



. (CIO) 



lL(i min ,zo) n f (L') 

Finally, inserting the above value of A coarse back into the expression for £ C oarse results in 

^coarse = A (5j , (Cll) 

i 

which is Adata- Thus, maximizing the measured cluster overdensity will give the correct richness 
and location of the cluster. 



C.2. Fine Grained C 

If the virtual bins are chosen to be sufficiently small that no bin contains more that one galaxy, 
then the calculation of C can be significantly simplified because there are only two probabilities 
that need to be computed. The probability of the empty bins 

^empty = e<*>** (C12) 

and the probability of the filled bins 

Pmicd = n J model e-<o dol . (C13) 
The sum of the log of the probabilities is then 

£fine = ^2 In -Pempty + ln jPfilled 

empty filled 

= - "model - "model + £ ln n model ( C14 ) 

empty filled filled 

By definition summing over all the empty bins and all the filled bins is the same as summing 
over all the bins. Thus, the first two terms in equation ( |C14 ) are just the total number of galaxies 
predicted by the model 



Yl n modei = / / n modci (d',l';z c )dQ'dl' 



all bins 



N f + AN C (C15) 
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where l m in is the apparent luminosity limit of the survey. Nf and AN C are the total number of 
field and cluster galaxies one expects to see inside 6> max ; they can be computed by integrating nj 
and n c 

/oo 
n f (l')dl' , 
-min 

N c {z c ) = / / n c (0',l';z c )dn'dl' 

Jlmin JO 

/"'"max 

= $ c (L(Z min ,z c )) / Z c (r')2irr'dr' (C16) 
Jo 

where $ C (L) = J£° <p c (L')dL'. 

Because we retain complete freedom to locate the bins wherever we like, we can center all the 



filled bins on the galaxies, in which case the third term in equation (C14) becomes 



n model 

^$>[r4 + A<] , (C17) 

filled i 

and the sum is now carried out over all the galaxies instead of all the filled bins. Combining these 
results we can now write the likelihood in terms that are readily computable from the model and 
the galaxy catalog 

£finc = -JV> -AN C + Y, ln[n) + An*] . (C18) 

i 

The simplest way to find the richness is to take the derivative of the above equation with 
respect to A and setting it equal to zero yields 



N, 



Si 



c 



E i T A * . ( C19 ) 



where we have substituted Si for n l c /n l j-. Unfortunately, it is not possible to solve for A nne explicitly, 
but a numerical solution can be found by standard methods. The resulting value of A nne is then 
inserted into the following equation for £ nne to obtain the maximized value of the likelihood 

£fine = -AfineA^ + £ ln[l + A flne ^] . (C20) 

i 

[Note: in the above expression terms that are independent of A nne and do not contribute any 
additional information to £ nne have been dropped.] 

Finally, it is worth mentioning again that while A coarse can be obtained directly from C coajTse , 
Afine can only be found by numerically finding the zero point of equation (C19). Furthermore, 
equation (C19) does not lend itself to standard derivative based solvers (e.g., Newton- Raphson) 
that produce accurate solutions in only a few iterations. Fortunately, the solution can usually be 
bracketed in the range < A nne < 1000, thus obtaining a solution with an accuracy AA ~ 1 takes 
log 2 (1000/1) = 10 iterations using a bisection method. 
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D. Data/Model Consistency Checks 

The first consistency check can be made with nf(l), or its cumulative probability 
Nf(l) = ff 30 n,f(l')dl', which can be fit directly to the galaxy catalog. In the case where a simulated 
catalog is used the luminosity distribution is taken from the underlying field luminosity function 
4>f(L) and it is possible to compute Nj(l) for each Hubble type directly 

N f (l) = 1^ $ f (L(l,z,»))d(z) 2 ^ldz , (Dl) 

where $/(/) = <t> f (l')dl' . 

The second consistency check can be made with the distribution in redshift (if estimated 
redshifts exist), where the number of galaxies at each redshift should satisfy the following 
convolution with the mean estimated error 
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Fig. 1. — Plots of the "2D" (i.e., a z — > oo) likelihood and richness as a function of redshift as 
computed from the coarse and fine matched filters. The input cluster has z c = 0.35 and a richness 
A = 100 (corresponding approximately to Abell richness class 1). 



800 




0.00 0.10 0.20 0.30 0.40 0.50 0.60 

redshift 

Fig. 2.— Plots of the "2±D" (a z ~ 0.05) likelihood and richness as a function of redshift as 
computed from the coarse and fine matched filters. The input cluster has z c = 0.35 and a richness 
A = 100. Although Afi ne is not peaked like A coarse this difference does not diminish the accuracy 
since the location of z c is determined entirely by £fi ne , which is sharply peaked. 
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a z ~ 0.05 z 
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Fig. 3. — Distribution of £ C oarse values. The range has been chosen so the distributions appear 
similar. The left dashed line denotes the value of the peak of the distribution. The right dashed 
line shows the value of £ cu t that was used. The horizontal solid line shows the FWHM around the 
peak. The significance levels of £ cut from top to bottom are approximately 8-a, 5-a and 5-cr. 
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Fig. 4. — Angular positions of all the simulated cluster galaxies with no background galaxies. 
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Fig. 5. — Angular positions of all the simulated cluster and field galaxies. 
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Fig. 6. — Richness and redshift of each of the simulated clusters. In this view the cluster parameters 
match the same 8 by 9 grid used in Figure ||. 
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Fig. 7. — Contour plot of £ C oarse computed from simulated data with photometric redshift errors 
a z ~ 0.05. Contour levels begin at £ C oarse = 40 and increase by a factor of 3 with each subsequent 
level. 
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Fig. 8. — Richness and redshift of each of the input clusters (boxes) with the short lines indicating 
the corresponding values determined from the AMF fine filter. The long curved line indicates the 
approximate detection limit and is computed by inserting £ cut = 40 into equation (12) and solving 
for A as a function z r - 
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Fig. 9. — Richness and redshift of each of the input clusters (boxes) with the short lines indicating 
the corresponding values determined from the AMF fine filter. The long curved line indicates the 
approximate detection limit and is computed by inserting C cn t = 100 into equation (12) and solving 
for A as a function z r . 
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Fig. 10. — Richness and redshift of each of the input clusters (boxes) with the short lines 
indicating the corresponding values determined from the fine filter. The long curved line indicates 
the approximate detection limit and is computed by inserting C cn t = 300 into equation (12) and 
solving for A function z r . 



N 



0.15 
0.10 

0.05 

0.00 

-0.05 

-0.10 

"8:1? 

0.10 
0.05 



^ 0.00 
N W -0.05 
-0.10 

-8:1? 

0.10 
0.05 
0.00 
-0.05 
-0.10 



-0.15 



0.00 



a ~ 0.05 



X 



X 



x~ 



X 



X 



X 
X 



i i i i i i i 



i i | i i i i i i i 

a z ~0.15 



X 



"1" 



X 



X~ 



X 



X 



X 



X 



X 



X 

X 



-X— 
X 



X 



I I I I I I I 



X 

_x 

X 



X 



X 



X 



X 



X- 

X 



I I I I I I I I I 

C z ->°o 



X 
X 



X 



X 



0.10 



0.20 



0.30 

Z 



0.40 



0.50 



0.60 



Fig. 11. — Redshift error as a function of redshift for all the detected clusters. Dashed lines 
indicate the \-a error range computed from these data, which from top to bottom are 0.014, 0.025 
and 0.023. 
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Fig. 12. — Fractional richness error as a function of redshift for all the detected clusters. Dashed 
lines indicate the 1-a range error computed from these data, which from top to bottom are 0.13, 
0.12 and 0.14. 



